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CHAPTER  I 
INTRODUCTION 


Juglans  nigra,  American  black  walnut,  is  the  most  valuable 
tree  species  in  Kansas.  Kansas  is  one  of  the  top  states  in  the 
nation  for  total  walnut  growing  stock  volume  and  growth  potential. 
It  is  the  second  ranking  state  in  domestic  walnut  veneer  log 
production. 

Black  walnut  is  common  in  the  eastern  one -third  of  the 
United  States  (Fig.  1).  The  range  extends  northward  to  southern 
Minnesota,  Wisconsin,  and  Michigan  including  the  southern  half 
of  New  York  to  the  Atlantic  coast.  It  reaches  southward  to 
northern  Florida  and  westward  to  eastern  Texas.  The  Mississippi 
Valley  and  Delta  regions  are  excluded,  but  it  does  well  in  the 
Appalachians  and  the  Midwest. 

The  natural  growing  range  of  black  walnut  extends  over  most 
of  the  eastern  one-third  of  Kansas.  This  area  contains  851  of 
the  1,564,000  acres  of  natural  woodlands  in  the  state,  and  receives 
the  most  rainfall.  In  the  western  one-half,  the  black  walnut  is 
generally  restricted  to  the  streams,  river  terraces,  and  drainage 
areas  where  moisture  is  sufficient  for  growth  (stricler  1973) . 


Figure  1.  Natural  Range  of  American  black  walnut  Juglans  nigra. 
(Reprinted  from  U.S.D.A.  Silvics  of  Forest  Trees  in 
United  States.  Handbook  No.  271.) 


Principally  the  stands  are  elm-ash-cottonwood  (bottomland 
type),  which  accounts  for  about  one-half  of  the  acreage.  The 
upland  forest  type  is  oak-hickory  and  is  usually  younger  than  the 
other  major  type. 

Kansas  is  estimated  to  have  about  60  million  cubic  feet  of 
black  walnut.  About  120  million  board  feet  are  sawtimber  logs 
and  62%  of  these  are  of  grades  1  and  2.  Black  walnut  has  been 
one  of  the  heaviest  cut  species  in  Kansas,  leaving  about  40%  as 
poletimber.  Poletimber  stands  provide  an  excellent  opportunity 
for  stand  improvement  techniques. 

Interest  in  planting  black  walnut  is  on  the  rise,  since  the 
value  of  the  crop  is  increasing  and  establishment  of  successful 
plantations  are  feasible. 

Site  selection  is  an  important  criterion  in  growing  black 
walnut.  The  soil  has  to  be  of  good  quality  to  yield  maximum  growth 
and  quality  logs.  Characteristics  such  as  drainage,  sufficient 
fertility,  thick  surface  horizon  of  silt- loam  or  loam,  and  neutral 
pH  must  be  present  to  obtain  a  fast  growing,  high  quality  crop. 
Best  sites  can  be  found  on  floodplains  and  terraces  or  protected 
coves,  but  naturally,  black  walnut  can  be  found  on  a  wide  variety 
of  conditions.  The  poorer  sites  usually  do  not  yield  a  high 
quality  log.  These  would  include  upland  sites  with  high  clay 
content  in  the  B  horizon,  and  tops  of  ridges  and  slopes  with 
shallow  soils. 

The  average  farm  has  many  potential  sites  for  growing  black 


walnut  that  are  too  small  to  be  productive  for  crops.  These 
include  edges  of  streams  and  small  coves  that  could  yield  valuable 
logs,  and  provide  other  benefits  to  the  landowner.  Knowledge  of 
the  potential  of  these  sites  would  be  of  great  benefit  to  urge 
the  landowner  to  invest  his  time  and  money  into  the  production  of 
black  walnut. 
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CHAPTER  II 
LITERATURE  REVIEW 


With  increasing  demands  for  timber  products,  foresters  are 
searching  for  ways  to  increase  productivity.  To  achieve  this  they 
have  a  way  to  evaluate  site  quality.  The  direct  estiamte  of  site 
index  is  the  most  widely  used  methods.  And  because  it  is  widely 
used,  more  information  has  been  developed  for  this  method  than 
for  any  other. 

With  the  appropriate  site  index  curves  for  a  particular 
area,  one  can  judge  (somewhat  accurately)  the  type  of  growth  to 
expect  from  a  particular  species  of  tree  for  that  type  of  site. 
An  example  of  direct  estimate  site  index  curves  is  a  graph  showing 
the  height- age  relationship  as  shown  in  Figure  2. 

After  a  site  index  graph  has  been  constructed  for  an  area 
and  species,  sites  can  be  evaluated  by  knowing  the  height  and 
age  of  dominant  trees.  The  two  values  are  found  on  the  graph 
and  the  point  of  intersection  lies  on  the  value  of  the  site 
index  curve.  Site  index  curves  are  usually  based  on  average 
stand  height  at  50  years  of  age. 

The  main  disadvantage  of  site  index  curves  are  that  they  do 
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Figure  2.  Site  Index  Curves  for  J.  nigra  by  Kellogg  (unpublished 
manuscript  covered  by  Brinkman  1966) 
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not  apply  to  all  sites.  Knowledge  of  the  site  index  for  one  area 
may  not  apply  to  another  (Jones  1969) .  Site  index  curves  must  be 
worked  out  for  each  species  and  area  to  be  evaluated. 

Other  factors  that  affect  direct  site  index  estimation  are 

these: 

1)  Stand  density  can  influence  tree  height. 

2)  Most  site  index  curves  are  based  upon  dominant- codominant 
trees,  while  curves  based  on  dominant  trees  alone  yield 

a  better  site  index. 

3)  A  very  large  sample  is  needed  to  avoid  error. 

4)  Tree  growth  varies  over  the  life  span  of  a  tree. 

5)  Genetically  superior  trees  are  utilized,  thus  not  giving 
a  true  representation  of  stand  potential. 

Making  a  site  index  curve  is  usually  time  consuming  and 
laborious.  In  the  standard  method  of  making  site  index  curves, 
trees  are  bored  and  a  ring  count  made  to  determine  age.  Total 
height  of  several  dominant -codominant  trees  are  plotted  over  average 
age  in  several  stands.  Site  index  curves  are  calculated  from  the 
regression  of  total  height  over  age.  Stem  analysis,  another  method 
of  direct  estimation,  is  accomplished  by  felling  the  tree,  taking 
a  section  out  of  the  tree  at  fixed  distances  up  the  trunk,  and  graphing 
a  series  of  height  over  age  coordinates,  thus  producing  the  tree 
curves  (the  stand  must  be  as  old  as  the  index  age  of  50  years)  for 
many  different  sites. 

An  alternative  approach  is  indirect  estimation  of  site  index, 
in  which  there  are:  vegetational ,  synecological  coordinates,  mensur- 
ational,  and  environmental  classifications. 


Vegetational  or  Plant  Indicator:  This  consists  of  classifying  the 
area  by  surveying  the  vegetation  and  usually  recording  climax  vegeta- 
tion. Much  work  has  been  done  in  this  area  but  with  little  useful 
knowledge  obtained  for  the  forester  (Jones  1969). 

The  Scandinavian  countries  worked  out  forest  habitat  as  compared 
with  climax-ground  vegetation,  soil  profile,  stand  structure,  and 
secondary  succession.  From  growth  studies,  yield  predictions  were 
made  (Jones  1969) . 

This  method  seemed  to  work  well  in  boreal  forests,  but  is  hard 
to  apply  in  an  area  such  as  the  United  States.  The  main  reasons 
for  difficulty  are: 

1)  The  geologic  material,  soil,  and  topography  are  sufficiently 
different  in  the  U.S.  that  they  should  be  used  as  evaluation 
criteria. 

2)  On  similar  type  of  soils,  the  tree  growth  alone  can  affect 
the  under  story  vegetation. 

3)  Trees  have  deeper  root  systems  and  are  affected  by  deeper 
soil  characteristics. 

4)  Stand  density  influences  the  climax  vegetation. 

5)  Key  plants  used  in  evaluation  are  not  visible  during  dormant 
seasons . 

Synecological  Coordinates :  Bakuzis  (1959)  developed  a  system  in 

Minnesota  using  vegetation.  First,  plants  were  rated  according  to 

their  needs  for  light,  moisture,  heat,  and  nutrients.  Then  the  site 

was  rated  by  the  presence  of  the  plants .  A  graph  was  developed  from 

the  data  showing  the  relationship  of  light,  moisture,  heat,  and 

nutrients  to  the  need  of  the  species  of  tree. 
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Mensurational  Methods :  Some  limited  success  has  been  achieved 
with  the  volume  growth  of  diameter,  and  aerial  photographs.  But 
this  is  used  when  other  methods  (Carmean  1975)  can  not  be  applied. 

Environmental  Aspect:  This  is  subdivided  into  the  holistic 
approach,  using  the  whole  environment  for  evaluation,  or  the 
factorial  approach  using  one  limiting  factor  as  the  key  to  a 
successful  apprasrial  of  the  site. 

Soil-site  evaluation  is  usually  considered  to  be  a  factorial 
approach  but  this  is  true  only  to  a  certain  point.  The  soil  is 
a  medium  through  which  outside  properties  can  act  to  affect  the 
quality  and  quantity  of  growing  space  and  conditions  for  roots. 
The  type,  structure,  and  properties  of  the  soil  regulate  the 
growth,  but  many  factors  influence  the  soil. 

Root  growth  is  primarily  determined  by  loil  moisture  and 
aeration  which  is  governed  by  the  soil  texture  and  percent  organic 
matter  (Coile  1952) .  Greater  depths  of  uniform  small  root  distri- 
bution is  achieved  as  the  coarseness  increases  in  the  soil.  In 
some  cases  stoniness  will  yield  better  growth. 

Some  studies  of  the  silt-plus-clay  content  of  the  A  and  B 
horizons  have  shown  that  there  is  decreased  growth  when  the  silt- 
plus-clay  increases  past  a  certain  percent  (Stoeckler  1948) . 
Plasticity  has  been  used  as  an  evaluation  guide,  with  about  the 
same  success  (Auten  1945b) .  This  is  not  surprising  since  plasticity 
is  influenced  by  the  clay  content  of  the  soil. 
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Other  characteristics  that  have  shown  good  correlations  have 
been  internal  drainage,  compactness  when  moist,  color,  depth  to 
mottles,  permeability,  and  thickness  of  A  horizon  (Coile  1952). 
The  texture  of  the  soil  influences  these  characteristics,  but 
texture  alone  gives  poor  correlations  to  the  site  index  since 
soils  of  the  same  textural  class  can  have  different  internal 
drainage,  aeration,  consistency,  and  structural  characteristics. 

Topography  has  to  be  observed  when  evaluation  a  site  because 
in  influences  moisture  availability  and  movement  (Auten  1945a) . 
Lower  sites  will  have  greater  potential  for  moisture  than  the  ridge 
tops  and  sides.  Slope  will  also  affect  the  depth  of  the  soil  to 
bedrock.  The  sides  of  ridges  usually  contain  shallow  soils  with 
greater  depths  at  the  bottom  of  the  slope  due  to  accumulation  by 
erosion  sediments.  Depth  to  bedrock  has  a  high  correlation  to 
site  index  (Auten  1937) .  Greater  productivity  has  been  found  with 
lesser  slopes  and  greater  depth  to  bedrock. 

Soil  fertility  and  acidity  has  been  found  to  be  the  least 
influential  of  the  soil  properties  (Ralston  1964) .  There  has 
been  no  one  nutrient  found  to  be  a  limiting  factor.  In  some  cases, 
increased  productivity  has  been  observed  with  potassium,  sodium, 
and  nitrogen.  But  usually  there  is  little  correlation  between 
fertility  and  site  index. 

Many  biological  factors  that  affect  tree  growth  have  to  be 
included  with  other  errors  in  evaluation  (Carmean  1975) .  These 
may  be  improper  stand  density,  competing  vegetation,  genetic  variations, 
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insects,  and  diseases.  Improper  stand  density  is  probably  the 
most  limiting  factor,  with  weeds  and  undesirable  species  robbing 
nutrients  and,  less  importantly,  water. 

Soil- site  evaluation  methods  have  yielded  high  correlations 
with  the  site  index.  Soil  characteristics  that  limit  soil  moisture 
and  aeration  correlate  the  best.  The  amount  of  clay  in  a  soil  is 
in  every  case  a  useful  evaluation  factor.  Percent  clay  will  influence 
the  plasticity,  cohesion,  porosity,  infiltration,  permeability,  and 
storage  capacity  of  the  soil. 

In  Kansas,  most  of  the  rain  comes  during  the  early  part  of  the 
growing  season,  with  little  stored  soil  water  accumulated  throughout 
the  winter.  Even  though  water  comes  when  the  plant  needs  it  the  most, 
there  is  not  enough  water  for  wasting  to  occur  and,  therefore,  good 
infiltration  and  storage  of  the  water  is  critical  for  good  growth. 

Depth  to  a  restricting  layer  such  as  a  claypan  or  fragipan 
or  a  dense  heavy  layer  governs  the  effective  thickness  of  the  soil 
(Ralston  1964).  As  the  soil  becomes  shallower,  the  yield  decreases. 
Recording  the  depth  and  strength  of  this  layer  could  easily  be 
done  in  the  field  and  would  provide  useful  information. 

Baker  and  Broadfoot  (1976)  developed  an  extensive  soil-site 
method  for  cottonwood  applicable  throughout  the  lower  Mississippi 
River  Valley.  They  evaluated  physical  condition,  available  moisture, 
nutrient  availability,  and  aeration  of  the  soil.  A  percentage  was 
given  to  what  each  contributed  to  the  maxumum  growth  at  the  age  of 
30.  These  four  major  factors  were  broken  down  into  minor  factors 
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categorized  by  their  level  of  contribution  to  the  major  factor. 
To  evaluate  a  site,  the  factors  were  measured  and  the  percentage 
of  maximum  growth  was  determined. 

The  study  reported  in  this  thesis  is  concerned  with  the  soil- 
site  approach.  Consideration  will  be  given  to  the  total  environment 
to  classify  the  material  that  is  collected  for  the  soil  characteristics, 


CHAPTER  III 
OBJECTIVES  OF  THE  STUDY 


The  objective  of  this  study  was  to  find  the  physical  soil 
characteristics  and  topographic  variables  that  affect  black  walnut 
tree  growth  in  northeast  Kansas.  The  usefulness  of  these  variables 
will  be  evaluated  by  devising  an  acceptable  mathematical  expression 
to  find  a  predicted  site  index  for  the  study  area.  This  predicted 
site  index  will  be  compared  to  the  observed  site  index  obtained 
from  the  conventional  height- over- age  site  index  curves  from  Kellogg. 

Variables  that  are  easy  to  measure  and  do  not  require 
highly  trained  personal  to  interpret  them  will  be  given  special 
consideration.  A  method  that  could  reasonably  estimate  the  site 
index  of  a  site  with  no  trees  present,  and  be  used  in  the  field 
by  non- specialized  personnel  would  be  of  great  value  to  the  field 
forester  and  potential  tree  farmer. 
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CHAPTER  IV 
METHODS  OF  INVESTIGATION 


Location  and  Description  of  Sites 

All  sites  were  in  the  northeastern  corner  of  the  state  of 
Kansas.  Counties  included  were  Pottawatomie,  Riley,  Geary, 
Marshall,  Nemaha,  Shawnee,  Jefferson,  Douglas,  and  Wyandotte  (Fig. 

3). 

Sites  were  selected  and  marked  by  W.  A.  Geyer  (Personal 
communication)  in  the  computation  of  his  direct  estimate  site 
index  curves  for  eastern  Kansas.  The  sites  were  selected  for  their 
presence  of  black  walnut  growing  in  a  natural  stand.  At  each  site 
dominant  trees  were  bored  with  an  increment  bore  at  stump  height 
(about  1/2  ft)  to  determine  their  age.  These  trees  were  selected 
to  avoid  open  grown  trees,  trees  with  broken  tops,  or  trees 
that  had  been  suppressed  by  other  trees.  They  were  then  marked 
with  tree  paint  and  numbered  for  later  identification.  Total 
height  of  each  tree  was  then  recorded  for  determination  of 
site  index  of  the  site  using  Kellogg 's  site  index  curves  for 
black  walnut  in  the  Midwest  (Fig.  2). 
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Figure  3.     Northeast  Kansas  showing  counties  used  and  location 
of  plots. 
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Of  the  sites  that  were  marked,  48  were  selected  for  the  soil- 
site  study.  A  base  age  of  50  years  is  the  ideal  age  for  trees 
because  that  is  the  base  age  for  the  Kellogg  site  index  curves, 
but  we  used  trees  from  age  11  to  111  due  to  time  and  location  of 
plots.  The  sites  were  selected  to  cover  a  range  of  site  index 
values.  Values  from  41  to  82  were  included  (Table  I). 

Field  Observations 

Each  site  was  visited  and  the  soil  profile  was  exposed.  A 
pit  was  dug,  12'  or  so  from  the  walnut  trees,  to  a  depth  of  3  to 
5'  depending  on  natural  barriers  such  as  rock.  If  a  natural  expo- 
sure was  present  it  was  utilized. 

The  horizons  were  then  identified  to  a  depth  of  5'  by  the 
standards  given  in  the  Soil  Survey  Manual  (1951  U.S.D.A.  Handbook 
No.  18).  If  horizons  were  not  present,  intervals  of  1'  were  used. 
Thickness  of  each  horizon  was  measured  to  the  nearest  inch.  Texture, 
structure,  boundary,  thickness,  and  rock  content  for  each  horizon 
or  1'  interval  was  recorded  (see  Fig.  4  for  form). 

The  texture  was  determined  by  the  "Feel  Method"  and  later  by 
laboratory  techniques  (Foth  et  al  1976) .  The  structure  and  boundary 
were  classed  according  to  standards  in  the  Soil  Survey  Manual. 
Rock  content  was  visually  observed  and  given  a  percentage  for  each 
horizon. 

The  depth  of  the  restricting  layer  was  then  determined  by 
measuring  how  deeply  the  small  roots  of  the  trees  penetrated,  and 
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Table  I.  Site  Indexes  and  Tree  Ages  for  48  sites, 


County 

Stop  No. 

Site  Index 
(Kellogg) 

Avg.  Tree  Age 

Pottawatomie 

2 

67 

21 

ii 

5 

59 

48 

M 

6 

52 

36 

H 

7 

65 

40 

ti 

8 

51 

34 

it 

9 

61 

32 

tt 

10 

49 

63 

it 

12 

45 

60 

ii 

13 

61 

90 

it 

15(1) 

48 

51 

ti 

15(2) 

48 

51 

M 

16 

40 

37 

M 

17 

60 

39 

it 

18 

58 

37 

it 

19 

63 

26 

n 

20 

48 

19 

H 

21 

49 

47 

ti 

22 

70 

60 

Riley 

1 

49 

48 

n 

3 

61 

29 

ii 

4 

62 

29 

18 


Riley 

5 

70 

26 

it 

6 

67 

'  27 

it 

7 

65 

101 

ii 

8 

70 

46 

Geary 

1 

43 

53 

H 

3 

41 

62 

ii 

4 

54 

53 

ii 

5 

63 

110 

n 

7 

48 

11 

Marshall 

2 

58 

58 

ii 

3 

75 

49 

n 

4 

50 

42 

ii 

6 

59 

62 

n 

7 

82 

57 

it 

8 

80 

40 

ii 

9 

62 

98 

Nemaha 

1 

62 

46 

Shawnee 

1 

75 

59 

Jefferson 

1 

79 

20 

ii 

5,6 

60 

78 

n 

7 

47 

32 

ii 

8 

60 

27 

Douglas 

8 

82 

34 

M 

9 

58 

54 

ii 

11 
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testing  the  strength  of  each  layer  by  the  penetration  of  a  knife 
blade.  The  restricting  layer  was  often  indicated  by  resistance 
to  shovel  penetration  while  digging  the  pit.  If  there  was  no 
restricting  layer  or  the  small  roots  were  throughout  the  profile, 
the  depth  was  recorded  as  5 ' . 

The  source  of  the  soil  parent  material  was  recorded  and  a 
visual  description  of  the  topography  was  made. 

Samples  of  every  horizon  or  foot  interval  were  collected  and 
labeled  for  laboratory  analysis.  A  lqt.  plastic  freezer  bag  was 
filled  with  soil  for  every  sample.  Naturally  occurring  clods  were 
then  taken  for  bulk  density  analysis  for  the  A  horizon  and  the 
restricting  layer,  if  present.  Clods  of  about  30  gr  each  were 
stored  in  plastic  bags  for  transportation  and  to  prevent  drying. 

The  pit  was  then  closed  and  restored  to  as  natural  a  condition 
as  possible. 

Laboratory  Procedures 

Bulk  density  was  determined  by  the  paraffin -clod  technique 
(Blake  1965) .  A  naturally  occurring  clod  of  about  30  gr  was  brushed 
to  remove  any  loose  material  and  examined  for  holes,  roots,  or  stones, 
If  any  were  found  the  clod  was  discarded.  A  thread  was  tied  around 
the  clod  and  the  clod  was  weighed. 

Paraffin  in  a  beaker  was  heated  in  a  water  bath  to  about  59*  C, 
enough  to  just  melt  the  paraffin.  The  clod  was  dipped  momentarily 
in  the  paraffin  several  times.  The  coating  was  checked  for  holes 
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and  air  bubbles.  Bubbles  were  broken  and  the  clod  redipped  if 
necewsary.  The  clod  was  then  weighed  with  the  paraffin  coating 
and  then  suspended  from  the  balance  and  weighed  in  water.  The 
water  content  was  determined  by  oven  drying  a  separate  sample  from 
the  same  horizon.  The  bulk  density  was  then  calculated  as  shown 
in  Appendix  A. 

The  texture  of  each  sample  was  determined  by  the  Bouyoucos 
Hydrometer  procedure  (Foth  et  al  1976) .  In  this  procedure  a  50  gr 
sample  of  soil  (100  gr  for  coarse  textures)  was  placed  in  a  Soil-Test 
baffled  stirring  cup.  The  cup  was  half  filled  with  distilled  water, 
5  ml  of  1  N  sodium  hexametaphosphate  was  added,  and  placed  on  a  mixer 
for  about  5  min  or  until  all  the  aggregates  were  broken  down. 

The  contents  of  the  cup  were  quantitatively  transferred  into 
a  Soil-Test  special  hydrometer  jar  and  was  filled  with  distilled 
water  to  the  lower  mark  of  1130  ml  (if  100  gr  was  used,  it  was 
filled  to  the  upper  mark  of  1205  ml) .  The  suspension  was  shaken 
horizontally  to  avoid  creating  circular  currents.  It  was  left 
undisturbed  for  20  sec,  the  ASTM  Soil  Committee  hydrometer  was 
inserted,  and  at  40  sec  the  first  reading  was  made.  The  suspension 
was  then  left  for  2  hrs  and  the  second  reading  made.  The  temperature 
was  recorded  before  the  first  reading  and  after  the  second.  A 
correction  factor  of  0.2  was  added  to  the  hydrometer  reading  for 
every  degree  above  or  below  67°  F.  for  each  reading. 

The  amount  of  sand  was  determined  by  subtracting  the  first 
corrected  reading  from  the  original  50  gr  sample  (the  hydrometer 
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reads  grams  of  soil  in  suspension) .  The  %   sand  was  obtained  by 
multiplying  the  amount  by  2.  The  second  corrected  reading  multi- 
plied by  2  yielded  the  I  clay.  Percent  silt  is  the  percentage 
remaining  from  1001.  Sample  calculations  are  in  Appendix  B. 

The  pH  of  every  sample  was  determined  by  using  a  Coleman  Portable 
pH  meter.  A  2  to  1  mixture  of  soil  and  distilled  water  was  placed 
in  a  paper  cup,  stirred  every  5  min  for  a  total  of  30  min.  The  pH 
was  then  immediately  read  after  the  last  stirring.  The  value  was 
rounded  to  the  nearest  tenth  of  a  pH  unit. 

Statistical  Analysis 

The  correlation  between  the  single  variables  and  the  site  index 
was  made  by  graphing  the  variables  against  the  site  index.  This 
was  done  for  the  variables  determined  in  the  field. 

Multiple  regression  using  the  SAS  76(5)  program  by  Barr, 
Goodnight,  Sail,  and  Helwig  (1976)  was  used  for  3  approaches.  First, 
all  the  quantitive  variables  were  entered  into  the  computer  and 
several  multiple  regressions  were  run  to  obtain  the  best  fit  of  the 
variables  to  site  index.  Stepwise  and  backwards  regressions  were 
used  to  find  the  best  variables. 

The  qualitative  variables,  SSOIL,  TOP,  nad  SURVEG,  could  not 
be  used  in  the  regression  model.  They  were  later  used  as  a  means 
for  sorting  the  sites  with  the  best  fitting  variables  to  obtain 
a  final  prediction  equation. 

Second,  the  sites  were  separated  by  SSOIL  and  a  regression 
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was  run  to  obtain  the  best  variables  for  each  (categories  3  (loess) 
and  4  (glacial)  were  entered  as  2  (colluvium-residuum)  for  this 
procedure) .  The  best  variables  were  used  with  the  SSOIL  to  obtain 
a  prediction  equation  for  each  source. 

Third,  the  easy  to  measure  variables  were  picked  and,  through 
trial  and  error  testing,  the  best  prediction  equation  was  formed. 
The  following  are  the  variable  names  generated  for  the  SAS  computer 
runs. 

Description  of  Variables  Tested 

DRL-  Depth  to  restricting  layer  in  inches  was  taken  in  the  field; 
the  procedure  was  described  under  field  observations. 

BDRL-  Bulk  density  of  the  restricting  layer  was  determined  in  the 
lab  for  2  samples.  The  value  tested  was  the  average  of  the  2  samples. 

TA-  Thickness  of  the  A  horizon  was  determined  in  the  field  to  the 
nearest  inch.  This  was  the  thickness  of  all  A  horizons  excluding 
and  A3  horizons  (these  were  considered  to  be  zones  of  transition) . 

BDA-  Bulk  density  of  the  Al  horizon  was  determined  by  laboratory 
procedures  for  2  samples.  Value  used  was  an  average  of  the  2  samples. 

SSOIL-  Source  of  soil  was  the  classification  of  the  soil  at  the 
site  as  1  of  4  types  as  listed  in  Table  II. 

TOP-  Topography  was  placed  in  one  of  the  categories  listed  in 
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Table  II.  Categories  used  in  qualitative  variables, 

SOURCE  OF  SOIL 

1  alluvium 

2  colluvium-residiuum 

3  loess 

4  glacial 


TOPOGRAPHY 


1  unproductive  dry 

2  gentle- rolling  terrain  (upland) 

3  upper  1/3,  top  of  slope,  or  ridge 

4  intermittent  stream  and  upland  cove 

5  flat  lst  terrace  (alluvial  bottoms) 

6  unproductive  wet 

7  stripmine 

8  middle  1/3  slope 

9  lower  1/3  slope 
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Table  II  using  the  description  of  the  site  in  the  field  notes. 

PCSLOPE-  Percent  slope  was  determined  in  the  field  by  the  Blume- 
Leiss  Altimeter. 

SURVEG-  Surface  vegetation  was  one  of  two  categories;  l=Duff 
composed  of  forest  litter  and  grown  vegetation,  and  2=Grass  with 
grass  predominantly  covering  the  ground. 

ESD-  Effective  soil  depth  (Steinbrenner  1965)  was  calculated  for 
each  profile,  effective  depth  being  the  depth  of  soil  minus  rock 
content.  First,  the  thickness  of  each  horizon  or  1  ft.  increment 
was  multiplied  by  the  percent  rock  of  that  horizon.  Then  the  inches 
of  rock  was  subtracted  from  the  thickness  of  that  horizon  or  1  ft. 
increment.  The  effective  soil  depth  was  the  total  of  new  values 
of  thickness  for  each  horizon  or  1  ft.  increment  for  the  total  5' 
profile . 

Ex.  Al  horizon  with  60%  rock  content  8"  thick 

8  x  .60  =  4.8    8  -  4.8  =  3.2  in.  for  ESD 

ESDRL-  Effective  soil  depth  to  restricting  layer  was  computed  to 
the  depth  of  the  restricting  layer  if  less  than  5 ' . 

WHCP-  Available  water-holding  capacity  (AWC)  was  given  in  inches 
of  water  per  inch  of  soil  for  every  profile.  This  was  accomplished 
by  building  a  chart  (Fig.  5)  to  estimate  the  AWC  for  each  textural 
class.  This  chart  is  a  combination  of  charts  by  O.W.  Bidwell  (1977), 
Salter  and  Williams  (1965),  and  Franzmeirer,  Whiteside,  and  Erickson 
(1960) . 
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Figure  5.  Textural  triangle  with  estimated  available  water  holding 
capacity  in  inches  of  water  per  inch  of  soil.  Adapted 
from  Bidwell  (1977) ,  Salter  and  Williams  (1965) ,  and 
Franzmeirer,  Whiteside,  and  Erickson  (1960). 
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None  of  these  charts  had  specific  values  for  each  textural  class. 
So  a  complete  chart  was  made  using  Franzmeirer,  Whiteside,  and  Erickson 
values  in  Michigan  (they  seemed  the  most  precise)  filled  in  by  values 
contained  in  the  other  charts.  Values  not  listed  in  the  3  charts 
were  estimated  by  studying  the  trends  of  all  3  charts. 

The  AWC  value  for  each  horizon  or  1  ft  increment  was  determined 
by  multiplying  the  new  thickness  value  for  each  horizon  or  1  ft 
increment  by  the  value  given  in  the  chart  for  the  texture  of  that 
horizon  or  1  ft  increment.  The  WHCP  is  the  total  of  the  AWC  values 
of  all  horizons  or  1  ft  increments  in  the  profile. 

WHCRL-  The  same  procedure  as  in  WHCP  was  followed  but  the  total 
was  stopped  at  the  depth  of  the  restricting  layer  if  less  than  5'. 

PHA1-  This  was  the  pH  value  to  the  nearest  tenth  pH  unit  for  the 
top  A  horizon  in  the  profile. 

PHB2-  This  was  the  pH  value  to  the  nearest  tenth  pH  unit  for  the 
B2  horizon.  If  the  horizon  was  subdivided  into  a  B21  and  B22,  the 
value  for  the  B22  was  used. 

PHC-  This  was  the  pH  value  to  the  nearest  tenth  pH  unit  for  the 
C  horizon  in  the  profile. 

LOWPH-  The  lowest  pH  value  in  each  profile. 

HIGHPH-  The  highest  pH  value  in  each  profile. 


SCB2-  The  %   silt  and  %   clay  were  added  together  for  the  B2  or 
B22  horizon  in  each  profile. 

SCA1-  The  %   silt  and  %   clay  were  added  together  for  the  top  A 
horizon  in  each  profile. 

PCCB2-  Percent  clay  of  the  B2  or  B22  horizon  in  each  profile. 

RATIO-  This  is  the  SCB2  value  for  1  profile  divided  by  the  TA 
value  for  that  profile. 

PCCA1-  Percent  clay  of  the  top  A  horizon  in  each  profile. 

ASP-  The  aspect  for  each  site  was  given  as  a  0-360°  compass 
reading.  This  is  the  angle  in  a  clockwise  direction  that  the  slope 
faced  from  north.  For  level  ground  the  value  was  0. 
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CHAPTER  V 
RESULTS  AND  DISCUSSION 


The  strongest  correlation  between  a  single  variable  and  site 
index  was  the  thickness  of  the  A  horizon  (Fig.  6).  The  correlation 
between  the  other  single  variables  and  site  index  was  low.  This 
was  determined  by  graphing  the  variables  against  site  index. 

The  backwards  elimination  procedure  proved  to  be  the  most 
useful.  It  started  with  all  the  variables  entered  and  used  them 
to  predict  site  index;  it  then  proceeded  to  eliminate  them  one  by- 
one  by  removing  the  least  significant  variable  with  the  highest 
PROB>F  (which  stands  for  the  lowest  alpha  level)  and  re-evaluating 
the  remaining  variables.  This  procedure  kept  eliminating  variables 
until  all  that  remained  were  significant  at  the  5%   level.  The  best 
fitting  variables  were  determined  by  myself  by  selecting  the  group 
with  the  highest  R- square  and  the  lowest  mean  square  which  is  the 
sum  of  the  difference  between  the  observed  values  and  the  predicted 
values  all  squared. 

The  stepwise  procedure  started  with  no  variables .  It  then  added 
the  most  significant  variable  with  the  lowest  PROB>F  and  determined 
the  R- square  for  that  one  variable.  It  then  looked  for  the  next 
variable  that  would  be  significant  given  the  f irstvariable .  If  a 
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Figure  6.  Graph  of  the  best  correlated  variable. 
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variable  was  entered  and  found  not  to  be  significant  by  5%  given 
the  previous  variables,  it  was  removed.  The  stepwise  regression 
proved  to  be  ineffective  because  it  worked  with  such  a  low  R-square 
(since  one  variable  explained  such  a  low  percent  of  the  model  state- 
ment) . 

The  variable,  bulk  density  of  the  restricting  layer  (BDRL) ,  had 
to  be  omitted  since  not  all  profiles  contained  a  restricting  layer 
and  because  the  PROB>F  was  high.  The  stepwise  and  backward  regres- 
sion were  run  with  this  variable  entered,  but  14  cards  were  thrown 
out  due  to  no  value. 

With  the  best  variables  obtained  (Table  III)  the  qualitative  varia- 
bles were  entered  in  as  classes  to  the  model  statement.  The  print 
out  contained  the  predicted  values  for  every  site  and  the  90%  confi- 
dence limits  for  each.  Four  data  cards  were  eliminated  from  the 
deck  for  they  contained  some  error  to  make  the  residual  deviate  more 
than  twice  the  standard  deviation. 

Two  were  for  trees  at  the  extreme  age  limits  for  acceptability, 
one  at  111  years  and  the  other  at  11.  At  these  ages  the  possibility 
for  the  site  index  curves  to  be  inaccurate  are  greater  than  if  the 
tree  was  close  to  the  base  age  of  50.  The  other  two  cards  were  for 
sites  that  had  very  high  site  index  values.  At  the  higher  values 
on  the  end  of  the  curve,  an  error  in  sampling  would  influence  the 
prediction  equation  to  deviate  the  curve  more  to  compensate  for  the 
error  than  if  it  was  located  along  the  middle  of  the  curve. 

With  the  removal  of  the  four  cards  the  stepwise  and  backward 
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multiple  regressions  were  rerun.  The  same  set  of  variables  in  the 
previous  analysis  plus  2  more  (BDA  nad  PHB2)  came  out  to  be  more 
significant.  The  mean  square  for  error  dropped  from  50.46  to  25.71 
and  the  R- square  increased  from  .64  to  .80  (comparison  of  Table  III 
to  Table  IV) .  Eleven  variables  out  of  the  20  tested  were  used  in 
this  procedure. 

With  the  11  quantitative  variables,  the  qualitative  variables 
were  examined.  This  was  done  by  trial  and  error.  All  possible 
combinations  of  the  three  variables  were  run  through  using  them 
to  sort  the  quantitative  variables.  If  the  PROB>F  was  large,  this 
indicated  that  the  variable  was  not  needed  given  the  other  variables, 

The  final  regression  is  shown  in  Table  V;  topography  was  the 
only  qualitative  variable  needed,  the  other  variables  came  from 
the  best  backward  regression.  The  R-square  of  .88  explains  that 
88%  of  the  variation  in  site  index  is  explained  by  the  model. 

The  prediction  equation  can  be  obtained  from  Table  VI.  The  topo- 
graphy of  each  site  must  be  rated  into  one  of  the  categories.  The 
equation  follows: 

SI=  Intercept  (11.77916362)  +  TOP  (times  the  value  given 
for  particular  catagory)  +  TA  (1.47156889)  +  BDA 
(12.26464736)  +  ESD  (0.47027600)  +  WHCRL  (4.19399357) 
+  WHCP  (-2.35752520)  +  PHA1  (10.52300802)  +  PHB2  (-3. 
23247513)  +  PHC  (6.82039023)  +  HIGHPH  (-9.88148137)  + 
PCCB2  (-0.48367093)  +  ESDRL  (-1.00747451) 
This  equation  had  a  90%  confidence  interval  of  about  17  units 
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(Table  VII).  The  highest  residual  was  -7.2,  with  20  of  the  44 
predictions  having  a  residual  of  less  than  2.0;  the  standard 
deviation  was  4.2. 

To  improve  the  equation,  another  approach  was  tried.  All 
the  sites  were  sorted  by  SSOIL  and  a  backwards  and  stepwise  was 
run  on  each  soil  source  yielding  a  set  of  variables  for  each  soil 
type.  Since  the  loess  category  had  only  one  site  and  the  glacial 
had  2  sites,  they  were  reclassified  into  the  colluvium-residuum 
category  putting  all  the  sites  in  2  categories. 

The  best  variables  for  each  soil  source  were  then  used  with 
the  source  of  soil  as  the  only  class  to  form  a  prediction  equation 
for  each  (Tables  VIII,  IX,  X,  XI,  XII,  XIII).  The  equations  are 
used  in  the  same  format  as  the  first  analysis;  the  predicted  value 
,is  the  total  of  the  sum  of  the  variables  multiplied  by  the  estimate 
of  each  with  the  intercept  added. 

The  regression  of  the  soil  variables  of  alluvial  soils  had  an 
R-square  of  .94  using  12  variables.  It  was  the  biggest  value  for 
any  procedure  tried.  The  90%  confidence  interval  was  about  14,  and 
the  mean  square  for  error  was  11.44.  The  largest  residual  was 
4.36  with  13  out  of  21  predictions  being  less  than  2.0  units  in 
deviation.  The  regression  of  the  best  variables  of  colluvium- 
residuum  soils  were  not  as  good,  the  R-square  being  .86  and  the  90% 
confidence  interval  of  about  19  units. 

An  ideal  method  would  involve  variables  that  could  easily  be 
measured  by  the  field  forester  with  no  laboratory  work  required. 
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With  this  in  mind  a  third  approach  was  tried,  with  the  variables 
picked  according  to  their  ease  of  measurement.  They  were:  SSOIL, 
TOP,  SURVEG,  ESDRL,  TA,  PCSLOPE,  PHA1,  and  ASP.  The  variables  were 
entered  in  the  model  using  the  qualitative  variables  as  classes. 

After  many  runs  with  different  combination  of  variables,  it 
was  found  that  TOP  was  used  as  the  class  and  TA  and  PCSLOPE  as  the 
other  variables  (Tables  XIV,  XV,  XVI) .  The  same  4  cards  as  in  the 
first  procedure  were  dropped  because  of  high  residuals.  This 
approach  had  a  R-square  of  .62  which  seems  low,  but  it  must  be 
remembered  that  only  3  variables  are  used  to  explain  62%  of  the 
model  statement  where  in  the  best  approach  using  soils  as  the  class 
and  getting  an  R-square  of  .94,  11  variables  had  to  be  observed. 

The  90%  confidence  interval  of  about  25  units  would  eliminate 
this  equation  from  being  used  in  the  field  for  this  spread  is  too 
great  to  be  used  on  black  walnut.  There  were  5  predictions  that 
had  a  residual  greater  than  10  units,  and  an  error  of  this  type 
could  cost  the  producer  greatly.  For  reliable  use  of  an  equation 
of  this  type,  you  should  have  a  901  confidence  interval  of  no  more 
than  about  10  units.  This  would  yield  estimated  site  indexes  that 
could  be  used  for  the  management  practices.  So  far  no  one  has 
developed  such  a  method  for  estimation.  This  study  using  the 
approach  of  classification  by  alluvial  soils  came  close  to  this  mark. 

As  found  in  this  study,  it  would  be  best  to  use  different 
approaches  for  different  types  of  sites.  If  the  site  was  of  alluvial 
source,  it  would  be  best  to  rate  it  using  the  prediction  equation 
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by  SSOIL  since  the  R-square  for  the  best  quantitative  variables 
classed  by  SSOIL  was  .94  for  alluvial  soils.  If  the  site  was 
anything  but  alluvial,  it  would  be  best  to  use  topography  to  class 
the  site  since  the  R-square  is  .88. 

Only  4  variables  came  out  significant  in  all  approaches  tested 
(easy  to  measure  approach  omitted  due  to  variables  being  hand 
picked)  they  were:  ESD,  WHCP,  PHC,  and  HIGHPH.  The  thickness  of 
A  horizon  (TA)  came  out  significant  in  all  approaches  but  for 
alluvial  soils.  This  is  not  surprising  since  most  of  the  alluvial 
soils  have  a  very  thick  A  horizon  (no  limiting  factor) .  Four 
variables  that  were  tested  were  never  significant,  they  were:  DRL, 
BDRL,  SURVEG,  AND  SCA1.  In  all  3  approaches  at  least  half  of  the 
variables  changed  significance  as  a  different  approach  was  tried. 

While  this  study  uses  sites  having  natural  growing  trees,  most 
studies  are  done  on  plantations  and  areas  where  one  can  obtain  many 
sites  on  the  same  soil  series.  This  study,  then,  is  influenced  by 
considered  as  errors.  I  tried  to  get  a  handle  on  this  natural 
variation  by  classifying  the  sites  by  topography  or  source  of  soil. 

One  of  the  problems  of  classification  by  topography  is  that  it 
was  extremely  difficult  to  decide  between  classifications  on  many 
occasions.  The  topography  did  not  always  fit  exactly  out  defi- 
nition or  in  many  cases  the  topography  would  be  changed  by  man's 
influence.  This  was  one  of  the  reasons  why  the  source  of  soil  was 
tried,  it  was  easier  to  class  the  site  as  alluvial  or  other.  It 
would  be  questionable  if  this  study  is  suitable  for  use  in  the  field. 
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Post  and  Curtis  (1970)  who  estimated  site  index  by  soils  and 
topography,  stated  that  their  method  would  not  be  useful  if  the 
standard  deviation  was  8.4  or  above.  The  highest  standard  deviation 
in  this  study  was  4.47,  but  even  with  this  range,  the  difference 
in  crop  value  for  black  walnut  would  be  great.  In  another  study 
Phillips  (1966) ,  for  yellow-poplar,  used  soil  and  topography  to 
estimate  site  index.  The  highest  R-square  achieved  was  .67  for 
7  variables. 

Many  soil-site  evaluations  (this  study  included)  are  made 
with  correlation  of  the  predicted  values  obtained  to  some  standard 
observed  value.  It  is  believed  that  this  method  is  not  a  reliable 
source  of  an  indicator  for  accuracy  of  the  prediction  (McQuilikin 
1976) .  A  better  method  would  involve  testing  the  acquired  prediction 
equation  on  independent  sites.  Cases  that  have  been  tested  this 
way,  have  found  to  be  completely  unreliable  for  a  site  index  predictor. 

There  are  many,  usually  interrelated,  variables  to  be  observed 
when  dealing  with  a  natural  situation  creating  many  problems  when 
attempting  an  estimation  of  the  results  of  these  variables  observing 
only  a  few  of  them.  It  may  be  possible  that  in  the  case  of  the  soil- 
site  approach,  the  soil  is  not  a  limiting  enough  factor  for  use  in 
obtaining  an  extremely  highly  significant  prediction  equation.  If 
this  is  gound  to  be  so ,  a  new  direction  must  be  followed  to  find 
the  answer. 


CHAPTER  VI 
SUMMARY  AND  CONCLUDING  REMARKS 


The  direct  estimate  of  site  index  that  this  study  used  as  a 
standard  for  comparision  might  be  one  of  the  largest  errors  in  the 
experiment.  These  site  index  curves  are  based  on  the  tree  age  50. 
Younger  stands  of  30  years  of  less  deviate  from  the  curves  because 
of  the  stand-age  soil-site  interactions  (Hannah  1971) .  Tree  growth 
on  loam  soils  might  be  very  rapid  at  first  and  taper  off  later, 
or  on  sandy  soils  tree  growth  would  start  slow  and  accelerate 
with  age. 

Stand  density  has  a  great  influence  on  tree  height.  Northeast 
Kansas  does  not  have  large  amounts  of  natural  black  walnut  occurring 
together.  Instead  the  walnut  trees  occur  in  small  patches  scattered 
throughout  the  woods.  Care  was  taken  in  selecting  the  sites,  but 
it  was  not  always  possible  to  avoid  low  stand  densities. 

Sampling  bias  is  a  major  possibility  with  harmonized  site  index 
curves.  Sampling  is  often  distorted  because  of  how  the  land  was 
originally  cut.  On  poor  sites,  old  trees  remain  and  the  good  sites 
contain  very  young  stands  (Beck  and  Trousdell  1973) .  This  is 
evident  in  northeast  Kansas. 
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The  second  error  in  harmonized  curves  is  the  assumption  that 
there  is  a  constant  curve  shape.  The  variance  to  site  index  can 
change  with  time  (or  age)  for  a  stand  (Lloyd  and  Hafley  1977) . 

Polymorphic  curves  developed  from  the  stem  analysis  approach 
could  yield  entirely  different  site  index  values  for  the  same 
site  (Beck  and  Trousdell  1973) .  The  problem  with  using  this  method 
on  black  walnut  is  that  the  trees  are  too  valuable  to  cut  down 
and  section  to  obtain  a  polymorphic  curve  for  the  site.  A  study 
of  this  type  could  be  carried  out  during  the  logging  of  the  black 
walnut  if  the  proper  cooperation  could  be  obtained  between  the 
loggers  and  the  researchers. 

Another  problem  encountered  in  this  study  was  the  procedure 
to  determine  bulk  density.  This  procedure  was  selected  because 
abundant  rocks  in  many  of  the  profiles  prevented  other  methods 
from  being  used.  The  clod  method  did  not  yield  the  same  values 
when  compared  to  values  obtained  by  other  methods.  The  values 
on  the  average  were  two  to  three -tenths  higher  than  what  was  expected, 
but  in  relation  to  each  other  the  values  seemed  to  be  consistent. 

The  higher  values  could  have  been  the  result  of  the  inter- 
clod  spaces  not  being  taken  into  account.  Another  reason  could 
be  drying  of  the  soils.  Soils  with  a  high  clay  content  shrink 
when  drying,  causing  the  air- dry  volume  to  be  less  than  the  field- 
moist  sample. 

This  study  is  a  very  positive  start  towards  soil-site  evaluation 
for  northeast  Kansas.  The  approach  of  separating  the  site  by  SSOIL 


and  obtaining  a  prediction  equation  for  each  type  shows  that  this 
could  be  done  for  topography  with  the  collection  of  more  sites. 
You  would  need  the  same  number  of  sites  as  the  number  of  variables 
tested  in  each  classification  for  backwards  and  stepwise  regresseion. 

With  the  re-evaluation  of  the  Kellogg  site  index  curves  by 
new  computer  techniques,  a  stronger  correlation  might  be  obtained. 
The  completion  of  Dr.  W.  A.  Geyer's  site  index  curves  for  northeast 
Kansas  could  improve  the  prediction  equation  further.  Both  these 
improvements  would  still  use  the  soil  information  collected  in  this 
study. 

This  study  showed  some  of  the  soil  characteristics  that  had 
a  strong  influence  on  the  site  index  value  of  the  site.  It  also 
eliminated  several  variables  from  being  tested  any  farther.  It 
was  a  strong  step  in  the  improvement  of  a  soil-site  evaluation 
technique  for  northeast  Kansas  and  with  further  work  a  field-working 
method  could  be  formed. 
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APPENDIX  A 
BULK  DENSITY  CALCULATIONS 
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1. 
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Moist  weight  -  dry  weight 

=  w 
dry  weight 


2.    Moist  weight  -  (  w  x  moist  weight  )  =  Adj.  dry  weight 


W 
3.    B.  D.  = 


Wa-Ww       Wa-W 
Pw     •    Pp 


B.D.=  bulk  density  of  soil 

W   =  Adj .  dry  weight  of  clod  before  paraffin  coating 

Wa  =  Adj .  dry  weight  of  clod  with  paraffin  coating 

Ww  =  Adj .  dry  weight  of  clod  with  paraffin  coating  in  water 

pw  =  Density  of  water  (1  gm/cc) 

p   =  Density  of  paraffin  (0.9  gm/cc) 
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APPENDIX  B 

SAMPLE  CALCULATION  OF 

MECHANICAL  ANALYSIS 
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SAND 


40  sec.  Corrected 

reading      Temp  «f      40  sec 

39.0        73.4         40.3 


(Temp  -  67)0.2  =  correction  factor 

(  73.4  -  67)  0.2  =  1.3 

(  50  -  corrected  40  sec.)  2  =  %  sand 

(  50  -  40.3)  2  =  19.41 


CLAY 


2  hr  Corrected 

reading      Temp         2  hr 


11.0         73.4         12.3 


2  hr  reading  x  2  =  %   clay 
12.3  x  2  =  24. 6* 


SILT 


100  -  %   clay  -  I  sand  =  %   silt 
100  -  24.6  -  19.4  =  56.0  I 


APPENDIX  C 
DATA 
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COUNTY         SJ5P  OBS    SI    DRL    BDRL    TA    BOA      SSOIL    TOP    PCSLOPE    SURVEC 

NO. 

POTT  2 

5 

6 

7 

8 

9 

10 

12 

13 

15(1) 

15(2) 

16 

17 

18 

19 

20 

21 

22 
RILEY  1 

3 

4 

5 

6 

7 

8 
GEARY  1 

3 

4 

5 

7 
MARSHALL       2 

3 

4 

6 

7 

8 

9 
NEMAHA  1 

SHAWNEE         1 
JEFF  1 

5,6 
7 
8 
DOUGLASS       8 
9 

11 
WYANDOTE       3 
4 


1  67 

34 

1.82 

16 

1.67 

4 

4 

0 

2 

2  59 

25 

1.70 

6 

1.50 

1 

5 

0 

2 

3  52 

13 

1.71 

13 

1.61 

2 

9 

10 

2 

4  65 

34 

1.90 

12 

1.56 

1 

4 

2 

2 

5  51 

17 

• 

10 

1.55 

2 

3 

20 

2 

6  61 

38 

1.7L 

15 

1.57 

1 

5 

0 

1 

7  49 

16 

1.71 

7 

1.51 

2 

3 

9 

2 

8  45 

29 

1.83 

4 

1.61 

2 

4 

15 

9  61 

60 

• 

12 

1.61 

1 

5 

0 

10  48 

18 

1.79 

5 

1.34 

2 

9 

100 

11  48 

29 

1.74 

6 

1.40 

2 

9 

100 

12  40 

11 

1.49 

11 

1.50 

2 

9 

80 

13  60 

27 

1.78 

13 

1.54 

2 

9 

40 

14  58 

34 

1.96 

18 

1.39 

2 

9 

15 

15  63 

32 

1.60 

10 

1.52 

1 

5 

0 

16  48 

12 

1.67 

5 

1.42 

2 

2 

5 

17  49 

23 

1.96 

10 

1.47 

2 

4 

5 

18  70 

60 

* 

16 

1.42 

1 

5 

0 

19  49 

24 

1.71 

15 

1.41 

2 

9 

2 

20  61 

15 

1.67 

8 

1.58 

2 

2 

15 

21  62 

20 

1.66 

6 

1.63 

2 

4 

15 

22  70 

19 

1.78 

12 

1.52 

2 

2 

10 

23  67 

29 

1.54 

13 

1.43 

1 

4 

0 

24  65 

34 

1.80 

7 

1.53 

2 

9 

100 

25  60 

32 

1.77 

8 

1.42 

1 

5 

0 

26  43 

17 

1.71 

10 

1.47 

2 

9 

10 

2 

27  41 

17 

1.67 

9 

1.46 

2 

9 

20 

28  54 

18 

1.83 

11 

1.45 

2 

9 

15 

29  63 

60 

• 

17 

1.47 

1 

5 

0 

30  48 

15 

1.70 

10 

1.41 

I 

5 

0 

31  58 

60 

• 

17 

1.47 

1 

5 

0 

2 

32  75 

60 

• 

21 

1.56 

1 

4 

0 

2 

33  50 

60 

« 

4 

1.80 

4 

8 

21 

2 

34  59 

60 

• 

11 

1.63 

1 

4 

0 

2 

35  82 

60 

• 

12 

1.39 

1 

5 

0 

1 

36  80 

60 

• 

23 

1.79 

1 

5 

0 

2 

37  62 

60 

« 

15 

1.63 

1 

5 

0 

38  62 

60 

* 

16 

1.77 

1 

5 

0 

39  75 

33 

• 

9 

1.64 

1 

4 

0 

40  79 

37 

2.01 
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A  soil-site  evaluation  was  made  from  forty-eight  plots  scattered 
throughout  northeast  Kansas.  These  plots  were  selected  for  their 
natural  occurring  black  walnut  growing  in  a  wooded  enviroment.  The 
trees  were  bored  and  measured  to  obtain  a  direct  estimate  of  site 
index.  The  site  index  curves  used  were  base  age  fifty,  from  Kellogg' s 

data  in  the  Midwest. 

At  every  site  the  profile  was  exposed  and  the  soil  horizons 
diagnosed.  Samples  were  collected  and  the  topography  of  the  site 

was  described. 

Physical  soil  properties  were  selected  as  quantitative  variables 
and  the  topography,  soil  source,  and  surface  vegetation  were  used 
as  qualitative  variables  for  sorting.  Three  approaches  were  tried: 
l./all  variables  on  all  sites  were  used  in  stepwise  and  backwards 
regression  to  obtain  a  prediction  equation,  2. /the  sites  were  sorted 
by  the  source  soil  and  a  prediction  equation  obtained  for  each,  and 
3. /the  easy  to  measure  variables  were  used  to  obtain  an  equation. 

The  predicted  values  were  compared  with  the  observed  values  for 
each  site  and  a  90S  confidence  interval  was  given  for  each.  The  best 
approach  was  sorting  by  alluvial  soils  with  a  R-square  of  .94.  With 
all  the  sites,  using  topography  as  the  class,  R-square  was  .88. 
Using  the  easy  to  measure  variables  proved  to  be  poorly  related  with 
a  R-square  of  .62;  this  approach  used  three  variables. 


